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INTRODUCTION 


The  detonation  of  nuclear  explosions  by  the  Soviet  Union  in  large  underground  cav¬ 
ities  under  either  a  Low-Yield  threshold  test  ban  treaty  (LYTTBT)  or  a  comprehensive 
test  ban  treaty  (CTBT)  would  constitute  the  greatest  challenge  to  U.S.  verification 
efforts.That  evasion  scenario  sets  the  limit  on  how  low  a  yield  can  be  verified  effec¬ 
tively.  The  recent  OTA  Report  Seismic  Verification  of  Nuclear  Testing  Treaties  found 
that  between  1  to  2  and  10  kilotons  (kt)  the  only  plausible  method  of  evasion  is  that  of 
nuclear  testing  in  large  cavities  in  salt  domes.  It  concluded  that  no  method  of  evading 
a  good  monitoring  network  is  credible  above  10  kt  and  that  several  evasion  scenarios, 
including  testing  in  cavities  in  bedded  salt  and  in  hard  rocks,  are  possible  below  1  to  2 
kt. 

Our  work  address  a  number  of  aspects  of  the  problem  of  clandestine  nuclear  testing 
in  large  cavities  in  salt  domes,  bedded  salt  and  hard  rock— what  types  of  evasion  sce¬ 
narios  are  plausable  based  on  geological  and  engineering  constraints,  which  ones  are 
likely  to  escape  detection  by  the  U.S.  and  which  ones  are  likely  to  be  identified.  We 
devote  particular  attention  to  the  Soviet  Union  and  to  the  critical  yield  regime  from  1  to 
10  kt,  where  scientific  research  over  the  next  few  years  seems  most  likely  to  have  a 
major  impact  on  the  verifiability  of  a  LYTTBT  and  on  what  threshold  can  be  verified 
effectively.  Several  important  aspects  of  decoupling  have  received  little  study  for 
more  than  20  years.  Since  then  a  great  deal  of  experience  has  been  obtained  by  the 
U.S.  and  several  European  countries  on  the  rheological  properties  of  salt  in  conjunction 
with  research  on  radioactive  waste  disposal  in  salt  deposits  and  by  industry  on  the  con¬ 
struction  and  stability  of  large  cavities  in  salt  for  petroleum  storage  and  waste  disposal. 
Experience  with  very  large  cavities  in  salt  domes  is  available  from  the  U.S.  Strategic 
Petroleum  Reserve  and  from  the  U.S.S.R.  and  U.S.  on  cavities  in  salt  for  gas  storage. 

Our  major  research  objectives  during  the  first  year  were: 

1)  Derive  revised  mb  values  for  Soviet  underground  nuclear  explosions  in  salt  and 
for  small  explosions  at  the  main  Soviet  test  sites;  determine  improved  magnitude-yield 
relationships,  especially  for  explosions  of  small  yield; 

2)  Study  numbers,  magnitudes  and  spectral  character  of  small  earthquakes  and  large 
chemical  explosions  in  and  near  areas  of  thick  salt  deposits  in  the  U.S.S.R.; 

3)  Evaluate  and  synthesize  data  from  engineering,  rock  mechanics  and  geological 
sources  on  the  properties  of  salt,  stability  of  large  underground  cavities  in  salt,  and  the 
use  of  cavities  in  salt  for  possible  clandestine  nuclear  testing,  especially  in  the  yield 
range  1  to  10  kt; 

4)  Start  work  on  maps  of  the  U.S.S.R.  showing  locations  suitable  for  possible  decou¬ 
pled  testing  of  explosions  of  various  yields  in  salt.  Identify  areas  that  have  known  salt 
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deposits,  those  that  are  not  known  to  but  conceivably  could  be  such  sites,  and  areas  that 
do  not  and  cannot  contain  salt  deposits  of  any  appreciable  thickness  (such  as  old  cra- 
tonic  regions  and  young  volcanic  areas). 

During  the  past  year  work  has  been  done  on  all  four  of  the  above  areas.  The  first 
topic  has  been  completed  and  a  paper  on  that  work  is  being  finished  for  journal  sub¬ 
mission.  Most  of  this  report  is  devoted  to  that  topic  and  to  the  implications  for  detect¬ 
ing  and  identifying  decoupled  nuclear  testing  in  the  Soviet  Union  in  the  Azgir  region 
and  in  adjacent  areas  to  the  north  of  the  Caspian  Sea  where  the  world’s  most  numerous 
salt  domes  are  found. 

There  has  been  a  long  debate  going  back  before  the  signing  of  the  TTBT  in  1974 
about  the  yields  of  Soviet  underground  nuclear  explosions  at  their  various  test  sites  and 
about  the  correct  methodology  for  calculating  explosion  yield  from  the  amplitudes  of 
various  seismic  waves.  The  establishment  of  a  threshold  in  terms  of  yield  by  the  TTBT 
put  a  premium  on  accurate  yield  determination.  After  1974  the  debate  in  the  United 
States  focused  upon  the  Shagan  River  portion  of  the  Soviet  testing  area  in  Central  Asia 
(eastern  Kazakhstan)  since  the  largest  underground  nuclear  tests  by  the  U.S.S.R.  after 
the  start  date  for  that  treaty  (April  1,  1976)  were  conducted  there.  (The  Shagan  River 
area  is  called  the  Balapan  testing  area  by  Soviet  workers.)  Until  three  years  ago  when 
Bocharov  et  al.  (1989)  published  yields  of  a  number  of  underground  explosions  at  the 
Central  Asian  test  site  for  the  period  1961  through  1972,  an  approximate  yield  had  only 
been  published  or  released  by  the  Soviet  Union  for  a  single  explosion  in  that  area.  The 
use  of  that  data  point  for  calibrating  yields  of  fully  contained  underground  nuclear 
explosions  was  difficult  since  that  event,  which  occurred  in  January  1965,  produced  a 
large  crater.  Its  mb  was  probably  smaller  than  that  of  a  fully  contained  explosion  of  the 
same  yield. 

Underground  nuclear  explosions  at  the  main  Soviet  test  sites  in  Central  Asia  and 
Novaya  Zemlya  generate  larger  P  waves  at  large  (teleseismic)  distances  that  U.S. 
explosions  of  the  same  yield  at  the  Nevada  Test  Site  (NTS).  The  difference  in  magni¬ 
tude,  often  called  the  mb  bias,  is  +  0.3  to  +0.4  magnitude  units  for  either  of  those  Soviet 
weapons  testing  sites  with  respect  to  NTS,  providing  the  explosions  are  in  the  same 
rock  type  (Office  of  Technology  Assessment,  1988).  In  practice,  most  U.S.  under¬ 
ground  explosions  have  been  conducted  in  softer,  less  competent  rocks  than  those  at 
Shagan  River,  leading  to  an  overall  bias  of  about  0.4  to  0.5  with  respect  to  NTS. 

In  1988  the  United  States  and  the  U.S.S.R.  conducted  nuclear  explosions  at  NTS 
and  Shagan  River,  respectively,  that  were  stipulated  by  bi-lateral  agreement  to  be  in  the 
yield  range  from  100  to  150  kt.  Each  country  permitted  the  other  to  make  on-site  mea¬ 
surements  of  yield  by  hydrodynamic  (Corrtex-type)  and  seismic  methods.  It  is  unfor¬ 
tunate  that  yield  estimates  and  other  pertinent  data  for  those  two  Joint  Verification 
Experiments  (JVEs),  as  well  as  yield  determinations  for  five  previous  nuclear  tests  at 
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each  site  that  were  communicated  to  the  other  country,  remain  classified  and  have  not 
been  announced  publically.  Nevertheless,  a  detailed  technical  report  on  the  two  JVEs 
in  the  New  York  Times  for  October  30,  1988  stated  that  the  American  and  Soviet  on¬ 
site  measurements  by  the  hydrodynamic  method  gave  yields  of  115  and  122  kt,  respec¬ 
tively  for  the  Soviet  JVE  of  September  14,  1988. 

In  previous  work  my  colleagues  and  I  determined  yields  for  Soviet  explosions  at 
Shagan  River  and  Novaya  Zemlya  by  using  calibration  curves  for  P  waves  that  were 
based  on  underground  explosions  of  announced  yield  for  which  their  P-wave  magni¬ 
tudes  were  corrected  for  mb  bias  between  the  actual  test  site  of  the  explosion  and  the 
Soviet  site  in  question.  Sykes  and  Ekstrom  (1989)  showed  that  the  reported  hydrody¬ 
namic  yields  of  1 15  and  122  kt  for  the  Soviet  JVE  were  consistent  with  those  mb- yield 
relationships  for  Shagan  River.  They  also  obtained  a  maximum  value  of  145  kt  for  the 
yield  of  that  explosion  using  surface  waves. 

In  the  last  few  years  much  progress  has  been  made  in  using  magnitudes  derived 
from  the  short-period  seismic  phase  Lg  for  precise  determinations  of  yield  of  explo¬ 
sions  at  Shagan  River  (Ringdahl  and  Fyen,  1988;  Ringdahl,  1989;  Ringdahl  and  Mar¬ 
shall,  1989;  Hansen  et  al.,  1990).  The  precision  of  a  single  determination  of  magnitude 
from  Lg  waves  is  often  comparable  to  that  obtained  after  averaging  values  of  mb  from 
50  to  100  worldwide  stations.  Magnitude-’  determined  from  Lg  waves,  m(Lg),  will 
probably  become  the  single  method  of  choice  in  the  future  as  more  data  become  avail¬ 
able  from  digitally-recording  stations  in  the  Soviet  Union  at  regional  distances.  Lg 
waves  studied  by  Ringdahl  and  his  co-workers  were  recorded  digitally  at  Norsar  and 
Graftenberg  at  distances  of  about  4000  km  from  Shagan  River.  Events  of  mb  less  than 
about  5.5  (Hansen  et  al.,  1990),  that  is  for  yields  smaller  than  about  40  kt,  typically 
have  a  small  signal-to-noise  ratio  at  those  large  distances.  It  appears  possible  to  use 
data  from  high-quality  stations  in  the  U.S.S.R.  to  determined  yield  from  Lg  waves  for 
tamped  (well-coupled)  explosions  as  small  as  1  kt  (Hansen  et  al.,  1990). 

Nevertheless,  accurate  determinations  of  yield  from  mb  continue  to  be  essential  for 
most  past  explosions  of  small  to  moderate  yield  at  Shagan  River  and  at  other  Soviet 
testing  areas,  explosions  detonated  in  the  1960s  before  digital  data  of  high  quality 
became  available  from  either  Norsar  or  Graftenberg,  Peaceful  Nuclear  Explostions 
(PNEs)  in  areas  of  the  U.S.S.R.  more  distant  from  those  arrays,  and  explosions  in 
regions  like  Azgir  where  the  propagation  of  Lg  is  poor  to  many  stations.  Yields  deter¬ 
mined  from  mb  become  a  valuable  check  on  calculations  from  Lg  and  can  be  used  to 
measure  the  precision  of  yield  determinations  by  the  two  methods.  I  show  later  that 
more  precise  yields  for  Shagan  River  explosions  can  be  obtained  by  suitably  averaging 
determinations  from  mb  and  m(Lg)  than  can  be  obtained  from  either  of  those  wave 
types  alone. 

A  major  problem  exists  in  converting  relative  yields  from  Lg  into  absolute  yields,  i. 
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e.  in  calibrating  a  given  test  site  in  an  absolute  sense.  While  yields  for  several  explo¬ 
sions  at  Shagan  River  prior  to  1973  have  recently  been  published  by  Soviet  workers 
(Bocharov  et  al.,  1989)  and  yields  of  five  other  Soviet  explosions  have  been  turned 
over  to  the  U.  S.  Government,  the  Soviet  JVE  of  1988  is  the  only  event  for  which  yield 
has  been  determined  independently  on-site  by  experts  from  outside  the  Soviet  Union. 
As  mentioned  above,  the  yields  measured  on-site  for  the  JVE  and  those  five  events 
have  not  been  formally  released  or  published  by  either  the  U.  S.  or  Soviet  govern¬ 
ments.  Thus,  calibration  of  absolute  yields  for  Shagan  River  from  Lg  alone  must  be 
based  on  either  the  single  data  point  for  the  Soviet  JVE  or  the  unverified  yields 
reported  by  Bocharov  et  al.  (1989). 

Body  wave-yield  relationships,  however,  have  the  advantage  that  data  from  many 
explosions  of  announced  yield  in  other  parts  of  the  world  can  “transported”  to  the 
Shagan  River  site  via  relative  bias  values  for  mb,  such  as  those  determined  from  atten¬ 
uation.  This  has  the  advantage  of  relying  on  more  than  a  single  mb-yield  pair  of  values 
in  constructing  a  magnitude-yield  relationship  for  explosions  at  Shagan  River. 
Ringdahl  and  Marshall  (1989)  calibrated  a  mb-yield  relationship  for  that  site  using  the 
mb  value  and  approximate  yield  (1 1 1  kt)  of  the  cratering  explosion  of  January  15, 1965 
and  then  making  a  correction  (addition)  to  its  mb  to  account  for  its  very  shallow  depth 
of  burial,  i.  e.  for  the  fact  that,  unlike  almost  other  Shagan  River  explosions,  it  was  not 
contained.  Since  Lg  measurement  were  not  available  for  that  event  from  Norsar  or 
Graftenberg,  Ringdahl  and  Marshall  (1989)  regressed  data  on  m(Lg)  with  those  on  mb 
for  many  contained  explosions  to  obtain  a  m(Lg)-yield  relationship  and  thence  to  com¬ 
puted  yields  of  other  events  from  that  expression.  Their  absolute  values  of  yield  are,  of 
course,  sensitive  to  the  uncertainties  in  the  values  of  mb,  yield  and  the  mb  correction 
term  for  the  1965  event. 

In  this  paper  I  compare  my  recalculated  values  of  mb  with  values  of  m(Lg)  reported 
by  Ringdahl  and  Fyen  (1988)  and  Ringdahl  (1989)  for  large  numbers  of  explosions  at 
Shagan  River.  The  difference  between  those  two  magnitudes  is  shown  to  vary  system¬ 
atically  from  southwest  to  northeast  across  that  testing  area.  That  variation  is  attrib¬ 
uted  to  systematic  differences  in  mb  (and  hence  in  mb  bias  with  respect  to  NTS),  with 
Lg  being  generated  uniformly  throughout  the  entire  Shagan  River  area.  Subdividing 
the  Shagan  River  area  into  three  parts,  corrections  to  mb  are  then  made  as  a  function  of 
location  on  the  test  site  to  form  a  corrected  body  wave  magnitude,  mb’  The  values  are 
normalized  to  the  southwestern  part  of  the  Shagan  River  area  where  mb  =  mb\ 

I  then  use  mb  values  for  explosions  of  announced  yield,  Y,  in  various  parts  of  the 
world  and  correct  them  to  the  southwestern  part  of  the  Shagan  River  testing  area  using 
the  bias  values  determined  from  the  relative  attenuation  of  P  waves  by  Der  et  al.  (1985) 
for  various  testing  areas.  A  regression  of  mb  upon  log  Y  is  obtained  for  those  data  but 
omitting  the  yields  of  Bocharov  et  al.  (1989).  Fortunately,  the  explosions  with  the 
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three  largest  yields  in  Bocharov  et  al.  were  located  in  the  three  different  subareas  of  the 
test  site  that  I  define.  When  their  mb  values  are  corrected  for  position  on  the  testing 
area,  the  corrected  values,  mb\  and  their  published  yields  fall  very  close  to  the  mb- 
yield  relationship  determined  without  recourse  to  those  data  points.  This  indicates  that 
the  Bocharov  et  al.  yields  must  be  nearly  correct  and  that  they  are  neither  seriously 
biased  nor  unrepresentative.  It  also  indicates  that  it  is,  in  fact,  mb  and  not  m(Lg)  that 
varies  systematically  with  position  on  the  test  site.  (Consistent  yields  are  not  obtained 
if  mbis  taken  to  be  constant  and  m(Lg)  varies  over  the  testing  area). 

Thus,  to  obtain  yield  estimates  for  Shagan  River  explosions  from  mb  that  are  more 
accurate  than  a  factor  of  1.4  ,  it  is  necessary  to  use  somewhat  different  mb-yield  cali¬ 
bration  curves  for  various  parts  of  that  test  site.  Values  of  m(Lg)  are  then  regressed 
against  mb’  so  that  a  new  magnitude  mb(Lg)  can  be  calculated  from  Lg.  In  that  way 
the  mb-yield  relationship  that  was  obtained  for  the  southwest  part  of  Shagan  River 
(Fig.  1)  can  be  used  to  obtain  very  accurate  yield  estimates  of  explosions  at  Shagan 
River  from  either  mb\  mb(Lg)  or  a  weighted  average  of  those  two  magnitudes.  An 
advantage  of  this  calibration  procedure  is  that  it  does  not  rely  on  the  magnitude  or  yield 
of  a  single  event  like  the  Soviet  JVE  or  the  1965  cratering  explosion;  it  also  permits  the 
yields  published  by  Bocharov  et  al.  (1989)  to  be  verified. 

The  yields  obtained  herein  by  an  equal  weighting  of  data  from  P  and  Lg  waves 
agree  very  closely,  i.  e.  within  about  15%,  with  those  calculated  by  Ringdahl  and  Mar¬ 
shall  (1989)  using  the  1965  cratering  event  for  calibration.  Seismology  appears  to  be 
clearly  capable  of  providing  yields  for  large  Shagan  River  explosions  that  are  at  least 
as  good  as  and  often  better  than  those  attributed  to  the  Corrtex  method  (Office  of  Tech¬ 
nology  Assessment,  1988).  For  small  explosions  at  that  site,  seismic  methods  have  a 
clear  advantage  for  yield  determination  compared  to  Corrtex  (Office  of  Technology 
Assessment,  1988).  For  other  less  well  calibrated  test  sites,  like  Novaya  Zemlya,  one 
to  a  few  Corrtex  measurements  can  be  used  to  calibrate  the  site  such  that  seismic 
waves  can  be  used  to  estimate  yields  more  accurately  for  past  and  possibly  future 
explosions  in  that  area. 


RESEARCH  ACCOMPLISHED 


Improved  Bodvwave  Magnitudes 

Revised  mb  values  were  recomputed  using  station  corrections  for  all  known  Soviet 
underground  nuclear  explosions  at  the  Shagan  River  testing  area  in  Central  Asia  and  in 
the  Azgir  region  to  the  north  of  the  Caspian  Sea  for  which  data  were  available  from  the 
International  Seismological  Center  (ISC)  and  the  United  States  Geological  Survey 
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(USGS)  for  the  period  January  1961  through  May  1988.  Station  corrections  for 
Shagan  River  events  were  derived  for  9  large  explosions  and  applied  to  all  undergroud 
events  for  which  data  were  available  for  that  testing  area.  Similar  previous  work  on 
magnitude  and  yield  determination  are  described  in  Sykes  and  Cifuentes  ( 1 984),  Sykes 
and  Ruggi  (1986,  1989),  Sykes  and  Davis  (1987)  and  Sykes  and  Ekstrom  (1989).  Sta¬ 
tions  used  in  the  recalculations  were  confined  to  the  distance  range  25°  to  95°. 

A  major  object  in  our  work  has  been  to  reduce  the  standard  error  of  the  mean  for  mb 
to  values  as  small  as  0.015  to  0.03  by  using  large  numbers  of  stations  (80  to  100  for 
larger  events)  and  by  applying  station  corrections.  Since  individual  stations  typically 
record  explosions  from  a  given  test  site  with  amplitudes  that  are  consistently  either 
higher  or  lower  than  the  mean  for  many  stations,  the  application  of  station  corrections 
considerably  reduces  the  standard  deviation  of  a  single  reading  and  avoids  biases 
related  to  the  inclusion  or  exclusion  of  individual  readings  from  one  event  to  another. 
It  is  particularly  important  to  use  station  corrections  since  data  from  the  many  stations 
of  countries  like  Canada  or  France  were  not  available  for  some  explosions  but  were 
available  for  many  others.  Similarly,  station  corrections  were  derived  for  the  seven 
largest  underground  explosions  at  Azgir,  where  the  testing  medium  is  salt,  and  applied 
to  other  events  at  Azgir,  including  the  PNE  explosions  of  1966  and  1968  for  which  the 
Soviet  Union  long  ago  announced  yields  of  1.1  and  25  kt  respectively.  Recomputed 
magnitudes  for  Azgir  and  Shagan  River  explosions,  their  standard  errors  of  the  mean, 
and  other  pertinent  data  are  listed  in  Tables  1  and  2. 

Small  Azsir,E*pl9.SiQns 

In  a  special  study  made  for  the  Azgir  explosion  of  1 . 1  kt  of  1966,  mb  was  measured 
at  1 6  stations  for  which  seismograms  were  openly  available,  giving  mb  =  4.524  +  .056. 
Although  small  in  amplitude,  P  waves  could  be  readily  identified  at  many  of  the  better 
WWSSN  stations  of  higher  gain  and  good  signal-to-noise  ratio.  The  ability  to  detect 
such  a  small  event  25  years  ago  using  analog  records  from  mainly  simple  (non-array) 
stations  reflects  the  high  coupling  of  a  tamped  underground  explosion  in  salt  and  the 
efficent  propagation  (high  Q)  for  P  waves  from  the  Azgir  area  to  stations  worldwide  A 
similar  study  for  the  1968  Azgir  PNE  event  gave  mb  =  5.529  ±  .027. 

Magnitude- Yield  Relationships  for  Shagan  River  and  Azgir 

The  mb  values  for  explosions  of  announced  yield,  Y,  in  various  parts  of  the  world  in 
either  hard  rock  or  below  the  water  table  were  normalized  to  the  southwestern  part  of 
the  Shagan  River  testing  area  and  to  Azgir  using  the  published  values  of  relative  atten¬ 
uation  and  mb  bias  of  Deretal.  (1985).  For  example,  a  bias  of  0.35 1  mb  units  between 
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the  southwest  part  of  the  Shagan  River  testing  area  and  NTS  was  used.  Values  of  mb  in 
Fig.  1  for  Soviet  explosions  are  from  Tables  1  and  2;  those  for  other  events  are  from 
Marshall  et  al.  (1979). 

An  aim  here  has  been  to  use  as  few  Soviet  announced  yields  as  possible  so  as  to  per¬ 
mit  the  yields  of  Bocharov  et  al.  (1989)  and  other  yields  recently  turned  over  by  the 
U.S.S.R  to  the  U.S.  to  be  verified.  The  resulting  mb-yield  regression  (dotted  line  in 
Fig.  1)  for  the  southwest  part  of  the  Shagan  River  area, 

mb  =  4.483  +  0.787  logY  (1) 

which  includes  yields  published  by  Bocharov  et  al.  (1989)  for  Shagan  River  ,  in  fact, 
differs  very  little  from  the  other  regression  in  Fig.  1  (solid  line)  that  does  not  include 
those  data  but  does  include  the  yields  of  PNE  explosions  at  Azgir  in  1 966  (1.1  kt)  and 
1968  (25  kt),  which  were  made  public  many  years  ago  as  part  of  an  exchange  of  data 
on  peaceful  uses  of  nuclear  explosions.  The  inclusion  or  exclusion  of  the  two  Azgir 
events  mainly  affects  the  slope  of  the  mb-log  Y  relation  and  not  the  computed  magni¬ 
tudes  for  yields  from  100  to  150  kt.  A  similar  relationship, 

mb  =  4.456  +  0.787  logY  (2) 

was  obtained  for  Azgir  using  a  bias  of  0.324  with  respect  to  NTS  (Der  et  al.,  1985).  It 
can  be  seen  that  the  two  explosions  in  salt  at  AzgiT  couple  as  well  as  those  in  hard  rock 
at  Shagan  River  as  do  the  three  in  granite  in  Nevada. 

Subdivision  of  Shagan  River  Testing  Area  for  Yield  Determination 

Magnitudes  determined  from  Lg  at  Norsar,  m(Lg),  and  our  recalculated  values  of 
mb  are  now  available  for  large  numbers  of  events  at  Shagan  River.  Fig.  2  shows  that 
the  difference  between  m(Lg)  and  mb  varies  systematically  across  the  Shagan  testing 
area  from  low  values  in  the  southwest  to  high  values  in  the  northeastern  part,  much  like 
that  found  previously  by  Ringdahl  and  Marshall  (1989).  In  constructing  Fig.  2  I  used 
the  locations  of  Marshall  et  al.  (1984)  as  well  as  Marshall,  written  communication,  for 
more  recent  events.  Using  m(Lg)  as  a  standard,  the  Shagan  River  test  site  is  divided 
into  three  subareas.  The  straight  line  boundaries  between  the  subareas  were  drawn 
such  as  to  minimize  the  number  of  misfitting  data  points  and  to  be  approximately  par¬ 
allel  to  structural  features  on  and  near  the  test  site.  A  new  magnitude  mb’  is  defined 
equal  to  mb  for  the  southwestern  area,  mb  +  0.145  for  the  northeastern  area,  and  mb  + 
0.067  for  the  transitional  zone  between  them.  Much  less  scatter  is  evident  for  the  sub- 
regions  and  for  the  entire  Shagan  River  area  (Figs.  3  and  4)  when  m(Lg)  is  compared 
with  mb’  rather  than  with  mb. 
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Yields  for  Shagan  River  Events  from  Lg  and  Body  waves 


A  new  magnitude, 


mb(Lg)  =  1 .090  m(Lg)  -0.440  (3) 

is  defined  from  the  regression  of  mb’  and  m(Lg)  that  is  shown  in  Fig.  4.  This  permits 
yields  for  Shagan  River  events  to  be  calculated  from  either  m(Lg)  or  mb’  using  (1). 
This  permits  Lg  observations  to  be  used  for  yield  determination  using  the  much  larger 
collection  of  mb  values  that  are  available  for  explosions  of  announced  yield  on  a 
worldwide  basis.  Its  use  also  avoids  using  mainly  unverified  Soviet  yields  when  calcu¬ 
lating  yield  directly  from  observations  of  Lg.  Yields  calculated  from  those  two  magni¬ 
tudes  are  much  more  consistent  for  the  same  explosion  than  those  calculated  from  Lg 
and  uncorrected  body  wave  magnitudes.  The  same  is  true  for  the  yields  of  the  three 
largest  explosions  published  by  Bocharov  et  al.  (1989)- 125,  140  and  165  kt,  which 
occurred  at  Shagan  River  in  1969  and  1972  in  each  of  the  subregions  defined  in  Fig.  2. 

The  quantities  mb\  mb(Lg)  and  the  average  of  those  two  magnitudes  were  each 
regressed  against  logY  using  the  published  yields  for  those  three  explosions  and  the 
average  yield,  119  kt,  determined  from  hydrodynamic  measurements  for  the  Soviet 
JVE  as  published  in  the  New  York  Times  (Sykes  and  Ekstrom,  1989).  The  slope  of 
each  magnitude-log  Y  relationship  was  held  constant  at  0.8  since  the  yields  cover  a 
very  narrow  range.  Each  of  those  magnitudes  had  the  following  standard  deviations: 
mb\  0.046;  mb(Lg),  0.035;  average  magnitude,  0.01 8  magnitude  units.  While  a 
slightly  better  weighted  average  of  the  two  computed  values  of  log  Y  could  be  calcu¬ 
lated,  only  a  simple  average  seemed  warranted  at  this  time.  Nevertheless,  it  is  clear  for 
these  four  events  that  the  precision  of  yield  estimates  can  be  improved  by  using  an 
average  based  on  the  magnitude  determined  from  Lg  and  the  corrected  magnitude  mb’ 
rather  than  relying  on  yield  estimates  based  on  either  of  those  magnitudes  alone. 

Close  agreement  is  found  between  the  yields  published  by  Bocharov  et  al.  (1989) 
and  those  calculated  from  the  relationship  in  Fig.  1  that  does  not  include  those  data  or 
any  other  recently  released  Soviet  data.  Thus,  the  published  of  yields  Bocharov  et  al. 
for  Shagan  River  must  be  equal  to  or  close  to  the  actual  yields.  The  inclusion  of  the 
four  yields  of  Bocharov  et  al.  for  Shagan  River  explosions,  in  fact,  leads  to  slightly 
larger  computed  yields  than  those  obtained  in  Fig.  1  without  those  data.  The  yield 
computations  based  on  mb’  in  Table  1  utilize  the  relationship  that  includes  those  data. 
Nevertheless,  they  do  not  differ  significantly  from  those  computed  with  the  other  rela¬ 
tionship  in  Fig.  1 .  The  close  agreement  between  yields  computed  from  mb’  and 
mb(Lg)  and  those  published  by  Bocharov  et  al.  also  argues  that  mb,  in  fact,  and  not 
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m(Lg)  varies  with  position  from  southwest  to  northeast  across  the  Shagan  River  testing 
area. 

Precision  and  Accuracy  of  Yield  Determinations  for  Shagan  River  Explosions 

The  values  of  log  Y  obtained  by  taking  a  simple  average  of  the  magnitudes  obtained 
from  Lg  and  body  waves  for  66  Shagan  River  explosions  of  Y  >  10  kt  (Table  1)  have 
standard  errors  of  the  mean  (SEM)  between  0.00  and  0.08  and  an  average  SEM  of  .029. 
These  translate  into  yield  determinations  based  on  those  two  wave  types  with  a  preci¬ 
sion  at  95%  confidence  (R  value)  of  a  factor  of  1 .0  to  1 .45  (mean  value  of  1 . 14)  (Fig.  5). 
Only  one  explosion  of  yield  greater  than  62  kt,  an  event  in  1969,  has  an  R  value  greater 
than  1.28.  The  decrease  in  precision,  i.  e.  increase  in  R  values,  for  events  smaller  than 
about  60  kt  probably  mainly  reflects  the  smaller  signal-to-noise  ratios  of  the  Lg  waves 
that  were  used  to  determine  m(Lg)  at  large  distances. 

These  uncertainties  in  yield  estimation,  of  course,  pertain  to  precision  and  not  nec¬ 
essarily  to  accuracy.  Nevertheless,  it  seems  unlikely  that  the  calibration  curves  used  in 
Fig.  1  are  more  uncertain  than  either  0.07  in  magnitude,  0.08  in  log  Y  or  a  factor  of  1 .2 
in  yield.  The  greatest  contribution  to  the  uncertainty  of  the  yields  determined  by  aver¬ 
aging  magnitudes  from  Lg  and  P  waves  for  explosions  with  yields  larger  than  about  60 
kt  comes  from  the  uncertainty  in  the  magnitude-yield  calibration  curve.  Nonetheless, 
yields  determined  in  this  way  for  explosions  larger  than  about  60  kt  appear  to  have  an 
accuracy  at  95%  confidence  that  is  comparable  to  and  perhaps  better  than  the  factor  of 
1 .3  attributed  to  the  Corrtex  method  for  Shagan  River  explosions  (Office  of  Technology 
Assessment,  1988).  The  published  yield  of  16  kt  (Bocharov  et  al.,  1989)  for  the  explo¬ 
sion  at  Shagan  River  of  February  10,  1972  agrees  well  with  the  yield  of  18  kt  (Table  1) 
obtained  from  mb\  This  indicates  that  the  calibration  relationship  and  yield  estimates 
for  that  test  site  near  20  kt  from  P  wave  magnitudes  probably  are  of  comparable  accu¬ 
racy  to  those  obtained  from  P  waves  alone  for  yields  between  100  and  150  kt.  At  the 
present  time  the  Corrtex  method  could  not  be  used  with  confidence  to  determine  yield 
for  explosions  smaller  than  several  tens  of  kilotons  (Office  of  Technology  Assessment, 
1988).  Thus,  seismic  methods  offer  a  clear  advantage  in  monitoring  compliance  under 
a  Low-Yield  Threshold  Test  Ban  Treaty.  The  availability  of  digital  recordings  of  Lg 
from  regional  stations  in  the  Soviet  Union,  such  as  stations  operated  in  conjunction  with 
the  IRIS  program,  should  result  in  greater  precision  in  determining  yield  for  tamped 
explosions  as  small  as  about  1  kt. 


Yields  Determined  from  Surface  Waves  for  Shagan  River  Explosions 

Yields  were  also  determined  in  Table  1  from  surface  waves  for  all  moderate  to  large 
nuclear  explosions  at  Shagan  River  for  which  digital  data  were  available  through  the 
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end  of  1986.  With  one  exception,  data  of  that  type  were  available  only  after  late  1978. 
Rayleigh  wave  (LR)  and  Love  wave  (LQ)  amplitudes  with  periods  in  the  narrow  band 
18-21  s  were  measured  on  vertical  and  transverse  component  digital  seismograms  for 
long-period  stations  in  Eurasia.  Only  data  for  pure  or  nearly  pure  continental  paths  were 
used  so  as  to  avoid  complications  resulting  from  surface  wave  propagation  across  con¬ 
tinental  margins.  The  method  of  analysis  for  computing  the  surface  wave  magnitude, 
Ms,  and  station  corrections  are  described  in  Sykes  and  Cifuentes  (1984)  and  Sykes  and 
Ekstrom  (1989). 

It  has  been  recognized  for  many  years  that  many  underground  explosions  at  Shagan 
River  trigger  the  release  of  natural  tectonic  stresses  in  the  vicinity  of  the  explosion  shot 
point.  Both  the  explosion  source  itself  and  the  release  of  those  stresses  contribute  to  LR 
and  Ms  whereas  only  the  latter  contributes  to  LQ.  Following  our  previous  work  on  yield 
estimation  from  Ms,  we  sought  to  make  a  correction  to  Ms  so  as  to  largely  remove  the 
effect  of  tectonic  release.  This  was  done  only  for  explosions  with  small  to  moderate  tec¬ 
tonic  release  compared  to  the  size  of  the  pure  explosion  itself,  i.  e.  for  LQ  /  LR  <  0.60, 
where  that  ratio  was  taken  as  an  average  over  all  stations  used  in  the  analysis.  A  full 
moment  tensor  solution  is  need  for  larger  values  of  LQ/LR. 

Yields  were  computed  from 

Ms +  B(LR/LQ)  =  2.16 +  0.97  log  Y  (4) 

using  the  measured  values  of  Ms,  where  the  second  term  is  a  linear  correction  to  Ms  for 
the  effect  of  small  to  moderate  tectonic  release  (Sykes  and  Cifuentes,  1984;  Sykes  and 
Ekstrom,  1989).  The  relationship  Ms  =  2.16  +  0.97  was  determined  by  Sykes  and  Cifu¬ 
entes  (1984)  using  explosions  of  announced  yield  in  various  parts  of  the  world  for  which 
either  LQ/LR  was  small  or  the  mechanism  of  the  tectonic  release  was  such  as  to  not  sig¬ 
nificantly  affect  Ms  as  averaged  over  a  range  to  azimuths  from  the  event.  Values  of  B 
were  obtained  for  each  explosion  at  Shagan  River  of  LQ/LR  <  0.60  for  which  Ms  was 
available  and  Y  had  been  determined  from  mb  and/or  mb(Lg).  The  average  value,  B  = 
0.45,  is  close  to  the  value  0.43  obtained  by  Sykes  and  Cifuentes  (1984)  for  a  much 
smaller  number  of  events  at  Shagan  River  for  which  moment  tensor  solutions  were 
available  at  the  time  in  the  literature.  That  value  is  in  accord  with  the  hypothesis  that 
tectonic  release  at  that  site  involves  predominantly  thrust  faulting  andjhence,  that  the 
correction  to  Ms  is  positive.  Also,  the  argument  can  be  turned  around  taking  B  =  0.43 
from  those  moment  tensor  events  to  derive  yields  from  surface  waves  that  differ  only 
slightly  from  those  in  Table  1 .  Those  yields  can  then  be  used  along  with  values  of  mb’ 
for  those  explosions  to  calculate  the  constant  in  the  mb-yield  relationship  for  Shagan 
River.  That  constant  is  nearly  identical  to  the  two  values  obtained  in  Fig.  1  using  a  bias 
value  of  0.35  mb  units  for  the  southwestern  part  of  the  Shagan  River  area  with  respect 
to  NTS. 

The  individual  values  of  log  Y  calculated  from  Ms  (Table  1),  however,  scatter  much 
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more  when  plotted  against  log  Y  determined  from  either  mb’  or  mb(Lg)  than  is  the  case 
when  the  latter  two  quantities  are  plotted  against  one  another.  Thus,  surface  wave  esti¬ 
mates  of  yield  are  less  precise  and,  hence,  are  not  used  in  this  paper  to  calculate  average 
yields.  They  can  be  used,  however,  as  a  less  accurate  additional  constraint  on  yield  if 
an  estimate  from  Lg  waves  is  not  available  or  the  average  yield  from  P  and  Lg  is  near 
or  above  the  150  kt  threshold  of  the  TTBT. 

Yields  of  Weapons  Deployed  on  Soviet  Strategic  Delivery  Systems  and  Tested  at 
Shagan  River 

In  November  and  December  1972  and  in  July  1973  the  Soviet  Union  detonated  three 
nuclear  explosions  at  Shagan  River  with  yields  between  140  and  200  kt.  Those  events 
occurred  well  before  either  the  negotiations  for  the  TTBT  in  July  1974  or  the  start  date 
for  that  treaty  in  1976.  The  U.S.S.R.  did  not  test  again  in  the  yield  range  above  120  kt 
at  Shagan  River  until  June  1979.  A  number  of  much  larger  explosions  with  yields  of 
500  kt  or  greater  were  detonated  at  the  two  remote  sites  in  Novaya  Zemlya  in  the  20 
months  between  the  signing  of  the  TTBT  and  its  start  date.  Those  three  events  at 
Shagan  River  were  the  only  nuclear  explosions  detonated  by  the  U.S.S.R.  at  any  of  their 
test  sites  in  the  yield  range  120  to  400  kt  from  1970  to  June  1979,  a  period  of  a  several 
years  both  before  and  after  the  signing  of  the  TTBT  (Sykes  and  Davis,  1987;  Sykes  and 
Ruggi,  1989). 

The  occurrence  of  three  explosions  at  Shagan  River  in  that  yield  range  within  9 
months  of  one  another  in  1972  and  1973  suggests  that  they  represent  the  full  yield  or 
nearly  full  yield  testing  of  important  and  distinct  Soviet  nuclear  weapons.  The  use  of 
precise  methods  of  yield  estimation  from  P  and  Lg  waves  permits  an  assessment  to  be 
made  of  the  yields  of  those  three  events  and  to  compare  them  with  the  published  yields 
for  the  two  events  in  1 972.  The  yields  determined  as  an  average  of  the  log  Y  values 
from  mb’  and  mb(Lg)  are  154, 138  and  193  kt  with  associated  R  values  (95%  confidence 
on  precision)  of  1.14,  1.22  and  1.22  kt  respectively  (Tables  1  and  3).  The  calculated 
yields  of  the  first  two  events  agree  closely  with  the  announced  yields  (Bocharov  et  al., 
1989)  of  165  and  140  kt  (Table  3).  The  explosion  of  July  23,  1973  of  Y  =  193  kt,  for 
which  the  yield  has  not  been  released  by  the  U.S.S.R.,  is  clearly  the  largest  underground 
explosion  at  Shagan  River. 

Those  three  explosions  occurred  a  few  years  before  the  initial  deployment  (Table  4) 
by  the  Soviet  Union  of  their  first  MIRVed  (multiple  independently  targetable  re-entry 
vehicles)  intermediate-range  ballistic  missiles  (IRBMs)  and  MIRVed  submarine 
launched  ballistic  missiles  (SLBMs)  according  to  Sykes  and  Davis  (1987).  The  throw- 
weights  of  the  SS-20,  (the  IRBM),  and  SS-N- 1 8  mod  1 ,  (the  SLBM),  and  the  number  of 
re-entry  vehicles  (three  per  missile)  are  consistent  with  warhead  yields  between  140  and 
200  kt  (Sykes  and  Davis,  1987).  SS-20  missiles  have  been  destroyed  under  the  terms 
of  the  Intermediate  Nuclear  Forces  treaty;  a  number  of  SS-N- 18  systems  are  still  oper- 
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ational  as  of  late  1991.  Thus,  one  or  more  of  the  three  largest  nuclear  tests  at  Shagan 
River  in  1972  and  1973  appears  to  have  been  for  the  weapons  to  be  carried  by  these  two 
important  additions  to  the  Soviet  arsenal  of  strategic  nuclear  weapons. 

Clustering  of  Yields  of  Shagan  River  Explosions 

The  newly  calculated  yields  in  Table  1  provide  strong  evidence  of  clustering  in  the 
distribution  of  yields  of  Soviet  tests  at  Shagan  River  from  the  start  of  the  TTBT  in  1976 
to  mid  1 988.  The  cumulative  number  of  events  at  Shagan  River  during  that  interval  with 
yields  less  than  Y  are  plotted  in  Fig.  6  as  a  function  of  Y.  The  yield  ranges  centered  on 
40,  88  and  138  kt,  where  the  slope  is  greatest  in  Fig.  6,  are  associated  wi'  the  most  fre¬ 
quent  testing  at  Shagan  River.  Additional  tests  with  yields  between  30  anti  :  00  kt  during 
that  period  occurred  at  Novaya  Zemlya  as  did  a  few  others  at  the  Degelen  Mountain  site 
in  Central  Asia.  Most  tests  smaller  than  about  30  kt  during  that  time  interval  took  place 
at  the  Degelen  Mountain  portion  of  the  Central  Asian  test  site  (Sykes  and  Ruggi,  1989). 

The  clustering  of  yields  near  138  kt  in  Fig.  6  undoubtedly  mainly  reflects  the  artificial 
cutoff  set  by  the  150  kt  limit  of  the  Threshold  Treaty.  Kidder  (1985)  shows  a  peak  in 
testing  by  the  United  States  from  1980  through  1984  at  yields  somewhat  smaller  than 
150  kt  that  is  quite  similar  to  the  Soviet  peak  near  138  kt.  He  states  that  the  “accumu¬ 
lation  of  [U.S.]  tests  near  150  kt  is  partly  the  result  of  testing,  at  reduced  yield,  strategic 
weapons  whose  yield  would  otherwise  exceed  the  TTBT  limit.”  Presumably,  a  similar 
logic  holds  for  at  least  some  Soviet  tests  near  138  kt  as  well. 


Compliance  with  the  Threshold  Test  Ban  Treaty 


As  yield  estimates  for  Shagan  River  explosions  have  become  more  precise,  the  num¬ 
ber  of  events  since  1976  with  a  yield  determined  from  either  P  waves,  Lg  or  their  aver¬ 
age  (Table  1)  that  exceeds  150  kt  has  become  smaller.  The  maximum  calculated  yields 
of  any  of  those  explosions  has  decreased  as  well.  For  example,  five  events  from  1976 
through  mid- 1988  have  average  yields  in  excess  of  1 50  kt:  157, 158, 162,  169  and  170 
kt  (Table  1).  Each  of  those  values  is  within  a  factor  of  1.13times  150  kt.  Only  the  event 
of  July  1973,  which  occurred  before  the  TTBT  became  effective,  has  a  yield  that  is 
greater  than  150  kt  at  a  high  level  of  confidence.  For  two  of  those  five  explosions  the 
average  yield  was  based  on  mh’  alone  since  Lg  magnitudes  were  not  available.  A  more 
common  situation  in  Table  1  is  for  one  of  the  two  yield  estimates  to  be  above  150  kt, 
another  below  and  the  average  somewhat  below  1 50  kt.  For  example,  the  event  of 
December  10,  1972  had  an  announced  yield  of  140  kt,  a  yield  from  mb  of  125  kt,  that 
from  mb(Lg)  of  153  kt,  and  an  average  calculated  yield  of  138  kt. 

Thus,  all  of  the  explosions  at  Shagan  River  from  the  start  date  of  the  TTBT  through 
mid- 1988  are  consistent  with  Soviet  compliance  with  that  treaty.  A  similar  conclusion 
was  reached  by  the  Office  of  Technology  Assessment  (1988).  Nevertheless,  the  new 
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estimates  are  still  not  accurate  enough  to  ascertain  a  10%  or  smaller  exceedance  of  the 
threshold  with  high  confidence.  For  example,  the  event  of  November  2,  1972  with  an 
announced  yield  of  165  kt  has  an  average  yield  of  154  kt  and  yields  from  P  and  Lg 
waves  of  144  and  165  kt.  If  any  events  since  1976  were  in  excess  of  the  threshold  by 
small  amounts,  the  most  likely  candidates  are  1 69  kt  on  August  1 8, 1979  and  1 70  kt  on 
October  27, 1984.  The  explosion  of  July  1973,  on  the  other  hand,  is  a  good  example  of 
present  capability  to  identify  a  past  or  future  event  with  a  yield  in  excess  of  about  190 
kt  as  being  significantly  larger  than  150  kt  at  a  high  level  of  confidence. 

Cavities  Formed  bv  Nuclear  Explosions  in  Thick  Salt  Deposits  that  Might  be  Usable  for 
Future  Clandestine  Nuclear  Testing 

The  U.S.S.R.  has  carried  a  number  of  nuclear  explosions  in  salt  near  Azgir  as  well 
as  conducting  several  other  PNEs  in  regions  known  to  contain  salt.  Cavities  produced 
by  nuclear  explosions  in  salt,  such  as  Salmon  and  Gnome  in  the  United  States,  have 
remained  standing  for  many  years  whereas  cavities  produced  in  other  rock  types  usually 
collapse  within  short  periods  of  time.  Many  investigators  working  on  decoupled 
nuclear  testing  that  might  be  conducted  under  either  a  CTBT  or  a  LYTTBT  have  paid 
considerable  attention  to  the  possible  use  of  cavities  produced  by  nuclear  explosions.  I 
will  show,  however,  that  monitoring  of  the  relative  few  areas  of  the  Soviet  Union  in 
which  cavities  of  that  type  could  exist  that  could  be  used  in  the  future  for  either  the  full 
or  nearly  full  decoupling  of  explosions  with  yields  larger  than  1  kt  is  a  tractable  prob¬ 
lem.  Instead,  more  attention  needs  to  be  devoted  to  the  feasibility  of  decoupled  testing 
in  the  yield  range  from  1  to  10  kt  in  large  cavities  produced  by  solution  mining.  Many 
more  potential  sites  are  available  in  the  U.S.S.R.  for  creating  large  cavities  by  solution 
mining  than  are  available  from  past  nuclear  explosions  in  salt. 

Table  2  lists  Soviet  underground  nuclear  explosions  that  were  detonated  through  1986 
either  in  or  in  the  general  vicinity  of  thick  salt  deposits.  As  mentioned  in  earlier  sec¬ 
tions,  yields  for  explosions  at  Azgir  were  determined  using  recalculated  mb  values  that 
included  station  corrections  derived  for  that  site.  Applying  those  station  corrections  to 
events  more  than  1 00  km  from  Azgir,  howe%  ',  did  not  reduce  the  standard  error  of  the 
mean  of  mb  values  for  those  events.  Hence,  nation  corrections  were  not  used  in  recom¬ 
puting  mb  for  the  other  events  in  Table  2.  The  yields  of  all  of  the  events  in  Table  2  were 
calculated  from  (2). 

Nine  of  the  explosions  at  Azgir  have  yields  between  25  and  110  kt.  All  of  those 
events  occurred  between  July  1 , 1 968  and  October  24, 1 979  in  the  area  of  numerous  salt 
domes  to  the  north  of  the  Caspian  Sea,  what  Soviet  workers  call  the  Pri-Caspian  region. 

I  will  now  make  estimates  of  the  sizes  of  decoupled  nuclear  explosions  that  could  be 
detonated  in  cavities  created  by  those  and  other  Soviet  nuclear  explosions  located  either 
in  or  in  the  general  vicinity  of  thick  salt  deposits.  The  conservative  assumption  will  be 
made  that  cavities,  in  fact,  remain  standing  for  all  of  those  events  and  that  they  could  be 
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used  for  conducting  small  decoupled  nuclear  tests  at  some  time  in  the  future.  First, 
information  about  the  depths  and  dimensions  of  the  cavities  created  by  the  U.S.  explo¬ 
sion  Salmon  and  two  explosions  at  Azgir  is  used  to  calculate  yields  of  nearly  decoupled 
(Sterling  conditions)  and  fully  decoupled  nuclear  explosions  that  could  be  conducted  in 
the  cavities  assumed  to  have  been  created  by  the  Soviet  events  in  Table  2. 

Salmon,  a  fully-tamped  explosion  of  5.3  kt,  was  detonated  in  a  salt  dome  in  the  state 
of  Mississippi  at  a  depth  of  828  m  in  1 964.  The  Sterling  nuclear  explosion  was  a  decou¬ 
pled  event  of  0.38  kt  detonated  in  the  Salmon  cavity  at  the  same  depth  in  1966  (Denny 
and  Goodman,  1990).  The  Salmon  cavity,  like  that  produced  by  the  Gnome  explosion 
in  bedded  salt  in  New  Mexico  and  that  produced  by  the  tamped  Soviet  explosions  of 
1966  and  1968  in  salt  domes,  was  not  perfectly  spherical  in  shape.  In  each  case  a  sig¬ 
nificant  amount  of  rubble  and  molten  salt  accumulated  at  the  bottom  of  the  initially 
nearly  spherical  cavity.  The  cavity  volume  for  the  Sterling  event  was  19,400  m3,  giving 
amean  radius  16.7  m  (Denny  and  Goodman,  1990).  The  1968  Azgir  event  of25ktwas 
detonated  at  a  depth  of  590  m  and  produced  a  cavity  with  a  volume  of  140,000  m3,  giv¬ 
ing  a  mean  radius  of  32.2  m  (Kedrovskiy,  1970).  The  1966  Azgir  explosion  of  1.1  kt 
was  detonated  at  a  depth  of  about  165  m  and  generated  a  cavity  with  a  volume  of  1 0,000 
m3  (mean  radius  of  13.4  m).  Since  that  cavity  filled  with  water  (Kedrovskiy,  1970),  it 
obviously  could  not  be  used  for  detonating  decoupled  explosions  unless  the  water  could 
be  removed.  The  Soviet  explosion  at  Orenburg  of  October  22, 1971  of  15  kt  in  bedded 
salt  was  used  to  produce  a  cavity  for  storing  gas  condensates  at  an  unspecified  depth. 
Its  volume  was  50,000  m3,  giving  a  mean  radius  of  22.9  m  (IzraeF  and  Grechushkina, 
1978). 

The  requirement  that  a  clandestine  nuclear  test  in  a  large  underground  cavity  in  salt 
not  produce  a  crater  or  other  disturbance  at  the  surface  and  that  it  be  fully  contained  so 
as  not  to  leak  radioactive  products  to  the  surface  set  limits  on  the  minimum  depth  of 
such  a  cavity.  For  a  very  weak  material  like  salt.  Latter  et  al.  (1961)  conclude  that  the 
pressure  on  the  cavity  wall  for  full  decoupling  must  be  less  than  or  equal  to  one  half  the 
overburden  pressure  (i.e.  the  vertical  stress,  pgh)  so  as  to  prevent  failure  in  tension  of 
the  surrounding  salt  material  and,  hence,  to  prevent  leakage  of  radioactive  gases  from 
the  cavity.  The  relationship  between  a  step  in  cavity  pressure,  P,  produced  by  an  explo¬ 
sion  of  yield,  Y,  in  a  cavity  of  volume,  V,  and  the  requirement  for  containment  that  P 
be  less  than  some  constant,  k,  times  the  vertical  stress  can  be  written 


P=  (y-l)Y/V<kpgh 


(5) 
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where  y  is  the  ratio  of  heat  capacities  taken  to  be  1.2  under  the  conditions  of  interest 
(Latter  et  al.,  1961),  p  is  the  average  density  of  the  material  from  the  surface  down  to 
depth,  h,  and  g  is  the  gravitational  acceleration  at  the  earth’s  surface.  The  average  den¬ 
sity  in  the  following  applications  is  taken  to  be  that  of  salt  at  the  Salmon  site,  2200  kg  / 

'l 

m  .  For  the  Latter  criterion  mentioned  above,  k  =  0.5. 

Taking  k  =  0.5  and  the  relevant  parameters  for  the  cavities  created  by  the  Salmon  and 
the  1966  and  1968  Azgir  explosions,  the  maximum  yields  of  fully  decoupled  explosions 
according  to  the  Latter  criterion  that  could  be  detonated  in  those  cavities  (after  convert¬ 
ing  the  energy  in  Joules,  J,  to  kilotons,  where  1  kt  =  4.184  x  1012  J)  are  0.21, 0.02  and 
1 .1  kt  respectively.  The  ratios  of  the  size  of  the  tamped  nuclear  explosion  that  was  used 
to  create  each  of  those  cavities  to  the  size  of  the  largest  fully  decoupled  event  that  could 
be  detonated  in  them  under  the  Latter  criterion  are  25,  55  and  23  for  the  Salmon,  1966 
and  1968  events. 

Sterling  was  almost,  but  not  fully  decoupled;  P  was  a  factor  of  1 .8  =  0.38  /  0.21  times 
larger  than  that  calculated  by  way  of  the  Latter  criterion.  The  value  k  =  0.90  calculated 
for  that  event  indicates  that  P  still  did  not  exceed  the  vertical  stress  (k  =  1).  Experience 
with  cavities  in  salt  used  for  gas  storage  indicates  that  leakage  of  the  storage  product 
may  occur  when  the  internal  pressure  in  the  cavity  exceeds  the  vertical  stress.  If  Ster¬ 
ling  conditions,  i.  e.  k  =  0.90,  apply  rather  than  the  Latter  criterion  to  the  1 966  and  1 968 
cavities  at  Azgir,  maximum  yields  for  explosions  that  could  be  detonated  in  those  cav¬ 
ities  are  0.04  and  2.0  kt.  The  ratio  of  the  size  of  the  tamped  nuclear  explosion  that  was 
used  to  create  a  cavity  to  the  size  of  the  largest  decoupled  event  that  could  be  detonated 
in  them  under  what  I  will  term  Sterling  conditions  is  13.9  =  5.3  /  0.38  for  the  Salmon  / 
Sterling  pair  and  27  and  12.5  for  the  1966  and  1968  explosions  at  Azgir. 

Thus,  much  larger  (12  to  55  times  larger)  tamped  explosions  are  needed  to  create 
cavities  than  the  maximum  sizes  of  either  fully  decoupled  or  nearly  fully  decoupled 
explosions  that  can  be  detonated  in  them.  This  is  important  since,  even  at  the  unclassi¬ 
fied  level,  the  United  States  knows  of  all  past  Soviet  underground  explosions  down  to  a 
small  yield  that  have  been  detonated  either  in  or  possibly  in  areas  of  thick  salt  deposits. 
The  yields  of  fully  decoupled  and  nearly  fully  decoupled  explosions  that  could  be  det¬ 
onated  in  cavities  that  may  remain  standing  from  those  events  is  at  least  an  order  of  mag¬ 
nitude  smaller  than  the  yields  of  the  explosions  that  generated  those  cavities. 

While  the  cavity  volume  produced  by  a  tamped  explosion  undoubtedly  varies  with 
the  depth  of  the  event,  it  should  be  appreciated  that  fully  decoupled  and  nearly  fully 
decoupled  nuclear  explosions  larger  than  1  kt  would  have  to  be  conducted  in  a  fairly  nar¬ 
row  range  of  depths.  An  air  filled  cavity  in  salt  is  likely  to  deform  significantly  at  depths 
greater  than  1000  m  whereas  an  evader  determined  not  to  be  caught  testing  at  such 
yields,  would  use  only  cavities  that  are  at  least  several  hundred  meters  deep  to  insure 
that  bomb-produced  products  did  not  escape  from  the  cavity.  There  is  a  tradeoff 
between  the  fact  that  a  shallower  tamped  explosion  produces  a  larger  cavity  than  a 


15 


deeper  one  in  the  same  material  and  the  fact  that  a  larger  cavity  at  a  shallower  depth  is 
needed  to  satisfy  inequality  (5)  for  a  decoupled  event  of  a  given  yield. 

Since  the  depths  of  most  of  the  explosions  in  Table  2  are  not  known,  two  estimates 
of  the  yields  of  fully  decoupled  and  nearly  fully  decoupled  nuclear  explosions  that  could 
be  detonated  in  cavities  created  by  those  events  have  been  made  in  the  last  column  of 
Table  2.  The  smaller  value  is  derived  by  dividing  the  computed  yield  of  each  event  in 
Table  2  by  the  yield  ratio  5.3  /  0.21  =  25,  where  5.3  is  the  yield  of  Salmon  and  0.22  is 
the  computed  yield  by  the  Latter  criterion  of  the  largest  fully  decoupled  yield  that  could 
be  tested  in  the  Salmon  cavity.  The  larger  yield  in  the  last  column  is  computed  by  divid¬ 
ing  the  calculated  yields  of  each  event  in  the  table  by  the  yield  ratio  for  the  Salmon  / 
Sterling  pair  =  5.3/ 0.38=  13.9.  Those  simple  calculations-- 1.0  and  1 .8  kt--for  the  cav¬ 
ity  produced  by  1968  Azgir  event  of  25  kt  are  very  close  to  the  values  obtained  above 
using  the  published  depth  and  volume  for  that  cavity— 1 .1  and  2.0  kt.  The  simple  calcu¬ 
lations  for  the  much  smaller  1966  event  overestimate  the  yields  calculated  above  from 
the  depth  and  volume  of  that  cavity. 

The  yields  of  the  various  events  in  Table  2  divided  by  13.9  (larger  numbers  in  last 
column)  are  taken  as  a  measure  of  the  yields  of  nearly  fully  decoupled  explosions  that 
could  be  detonated  in  the  cavities  produced  by  those  events.  This,  of  course,  assumes 
Salmon  /  Sterling  conditions.  Uncertainties  in  the  yields  of  the  events  that  produced 
those  cavities  and  their  depths  can  lead  to  uncertainties  in  the  estimates  of  decoupled 
yields  by  a  factor  of  1.5.  Our  main  concern  here,  however,  is  what  are  the  “opportuni¬ 
ties”  that  may  be  available  to  test  weapons  clandestinely  of  certain  approximate  sizes 
and  not  that  exact  size. 

Table  5  summarizes  the  yields  of  nearly  fully  decoupled  (Sterling  conditions)  nuclear 
explosions  that  could  by  detonated  in  the  cavities  of  events  in  Table  2  both  by  area  and 
yield  range.  Most  of  the  “opportunities”  for  such  testing  are  concentrated  in  the  area  to 
the  north  of  the  Caspian  Sea.  The  “oppontunities”  for  the  larger  tests,  up  to  a  maximum 
of  7  to  8  kt,  are  mostly  confined  to  Azgir  itself.  Possibilities  for  such  decoupled  testing 
are  more  restricted  in  the  area  of  bedded  salt  to  the  northwest  of  lake  Baikal.  Some  and 
perhaps  all  of  those  events  may  not  have  been  conducted  in  salt  but  in  other  rocks  in  the 
strategraphic  sequence  or  in  rocks  just  outside  the  area  of  salt  deposits.  Once  again,  this 
list  attempts  to  be  conservative  in  its  cataloguing  of  potential  areas  and  sites  for  decou¬ 
pled  testing  in  the  range  1  to  10  kt.  For  example,  the  single  Bukhara  event  (Tables  2  and 
5)  was  used  to  put  out  a  fire  in  an  oil  well  that  had  encountered  an  unanticipated  fault 
that  had  provided  pathways  for  escaping  petroleum.  It  was  also  detonated  at  such  a 
great  depth  (2.5  km)  that  closure  of  the  salt  cavity  undoubtedly  occurred  soon  after  it 
was  detonated.  Thus,  on  both  of  those  grounds  it  appears  to  be  a  poor  and  risky  candi¬ 
date  site  for  conducting  a  clandestine  decoupled  nuclear  explosion. 

Thus,  the  number  of  cavities  produced  by  past  nucler  explosions  that  potentially 
could  be  used  for  clandestine  testing  of  nearly  fully  decoupled  nuclear  explosions  under 
either  a  CTBT  or  a  LYTTBT  is  limited.  Those  sites  are  confined  to  a  few  areas  of  the 
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U.S.S.R.  The  question  of  using  those  cavities  instead  for  partially  decoupled  events  as 
described  by  Stevens  et  al.  (1991)  will  be  dealt  with  in  a  separate  paper  along  with  the 
feasibility  of  constructing  and  using  cavities  in  hard  rock  for  either  the  full  or  the  partial 
decoupling  of  explosions  with  yields  from  1  to  10  kt.  Possibilities  of  clandestine  testing 
in  large  cavities  created  by  solution  mining  also  will  be  dealt  with  in  a  separate  paper. 

CONCLUSIONS  AND  RECOMMENDATIONS 


Implications  for  Identifying  Small  Decoupled  Nuclear  Explosions  in  Salt 

For  purposes  of  appreciating  the  detection  capability  of  a  given  seismic  network,  it 
is  important  to  recognize,  using  data  from  Azgir,  that  a  fully-coupled  (tamped)  explo¬ 
sion  of  1  kt  in  salt  in  high-Q  areas  of  the  U.S.S.R.  has  an  mb  of  4.46;  fully  decoupled 
events  of  1  and  10  kt  have  mb’s  of  2.61  and  3.40  respectively  (assuming  a  decoupling 
factor  of  70).  These  magnitudes  are  higher  than  has  generally  been  thought  previously. 
For  example,  chemical  explosions  of  mb  <  2.6  in  high  Q  areas  containing  salt  need  not 
be  considered  in  monitoring  a  1  kt  threshold  treaty.  Most  areas  of  thick  salt  deposits  in 
the  U.S.S.R.  are  typified  by  high  Q  (efficient  transmission)  for  P  waves  and  low  natu¬ 
ral  seismic  activity. 

Past  nuclear  explosions  conducted  in  salt  by  the  Soviet  Union  for  which  cavities  may 
remain  standing  that  are  large  enough  for  the  full  decoupling  of  explosions  larger  than 
1  kt  are  concentrated  in  only  a  few  areas  of  that  country.  The  existence  of  all  cavities  of 
that  size  or  larger  that  have  been  created  by  past  nuclear  explosions  is  known  (Table  3) 
since  the  explosions  that  created  those  cavities  must  be  about  25  times  larger  in  yield 
than  the  size  of  the  fully  decoupled  event  that  can  be  detonated  in  them.  (That  ratio  is 
13.9  for  the  Salmon/Sterling  pair  where  the  yield  of  Sterling  was  about  1.8  times  larger 
than  that  of  a  fully  decoupled  explosion  in  the  Salmon  cavity). 

Hence,  the  monitoring  of  cavities  of  that  size  that  may  remain  standing  that  were 
created  by  past  nuclear  explosions  should  be  relatively  easy  under  a  Low- Yield  Thresh¬ 
old  Test  Ban  Treaty  with  a  threshold  of  1  kt,  providing  U.S.  stations  are  allowed  to  oper¬ 
ate  under  the  treaty  near  the  epicenters  of  those  past  explosions.  Probably  the  greatest 
difficulty  in  monitoring  such  a  LYTTBT  involves  cavities  created,  not  by  nuclear  explo¬ 
sions,  but  by  solution  mining  in  other  areas  of  thick  salt  deposits  of  the  U.S.S.R. 

Precise  Yield  Estimates  for  Explosions  at  Shagan  River 


Yield  determinations  for  Shagan  River  explosions  can  be  improved  by  subdividing 


that  testing  area  into  three  parts  with  different  mb-yield  relationships  and  by  combining 
those  results  with  yields  calculated  from  Lg.  The  calculated  yields  of  many  of  those 
explosions  since  1976  are  clustered  in  a  few  specific  ranges  of  yield. 
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Magnitude  mb  normalized  to  SW  area  of  Shagan  River  Testing  Area 
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Fig.  1  Calibration  curves  used  for  yield  estimation  as  normalized  to  southwestern 

portion  of  Shagan  River  testing  area.  Explosions  at  Degelen  and  Murzhik  were 
not  used  in  two  regressions.  Explosions  used  are  those  of  published  yield  in  either 
hard  rock  or  below  water  table. 
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Fig.  2  Lg  magnitude  minus  body  wave  magnitude  as  a  function  of  location  within 
Shagan  River  testing  area.  Testing  area  is  divided  into  three  parts  for  better 
estimates  of  yield  from  body  waves. 


Fig.  4  Lg  magnitude  as  a  function  of  revised  body  wave  magnitude,  mb'.  SD  =  standard 
deviation. 
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Table  2.  Underground  Nuclear  Explosions  in  or  Near  Thick  Salt  Deposits  of  U.S.S.R. 
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Table  3:  Largest  Shagan  River  Underground  Explosions 
prior  tojune  1979  of  Yield  >  125  kt 


Date 


y.mm 


02  Nov.  1972 

165 

154 

10  Dec.  1972 II 

140 

138 

23  July  1973 

— 

193 

Table  4:  Major  Soviet  Missile  Systems  with  Throw-weights 

Appropriate  to  above  Yields 


System  First  Deplovii 

SS-20  1977 

SS-N-18  mod  1  1978 


#RVs  Comments 

3  First  Soviet  MIRVed  IRBM 
3  First  Soviet  MIRVed  SLBM 
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Table  5.  Number  of  cavities  that  may  remain  standing  from  Soviet  Underground  Nuclear 
Explosions  form  1961  to  1986  within  and  near  salt  deposits  that  might  be  useful  for 
conducting  nearly  fully  decoupled  nuclear  explosions  with  yields  greater  than  0.5  kiloton. 
(Obtained  from  calculated  yield  of  known  explosion  at  site  divided  by  yield  ratio  for 
Salmon/Sterling  explosions  =  5.3/0.38=13.9.) 


Area  of  Known  Salt  Deposits 


1.  North  Caspian  Region 

a.  Azgir 

b.  Astrakhan 

c.  Other 


a.  Central  Siberian  Platform 
(northwest  of  Lake  Baikal) 

b.  Within  a  few  hundred 
kilometers  of  layered  salt 
deposits 


3.  Bukhara  II,  Central  Asia 
(explosion  used  to 
extinguish  fire  in  oil  well) 


Maximum  Nearly  Fully  Decoupled  Yield  (kt) 


Many  small  nuclear  explosions  at  nearly  the  same  location.  Cavities  conceivably 
could  be  interconnected  to  form  a  larger  cavity 
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